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Abstract. We study the overdamped dynamics of a Brownian particle in the double- 
well potential under the influence of an external periodic (AC) force with zero 
mean. We obtain a dependence of the jump rate on the frequency of the external 
force. The dependence shows a maximum at a certain driving frequency. We 
explain the phenomenon as a switching between different time scales of the system: 
interwell relaxation time (the mean residence time) and the intrawell relaxation time. 
Dependence of the resonant peak on the system parameters, namely the amplitude 
of the driving force A and the noise strength (temperature) D has been explored. 
We observe that the effect is well pronounced when A/D > 1 and if A/D ^> 1 the 
enhancement of the jump rate can be of the order of magnitude with respect to the 
Kramers rate. 
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1. Introduction and the background 

In the last decades a lot of interest has been devoted toward the problem of the behaviour 
of a nonlinear system under the combined influence of stochastic and time-periodic 
forces. A number of remarkable phenomena such as stochastic resonance (SR) and 
resonant activation has been discovered and extensively studied both experimentally 
and theoretically [1, 2, 3]. So far, a lot of experimental confirmations of this effect has 
been reported in different areas of physics, like optics, biophysics or condensed matter 
to name a few. In our opinion, the frequency dependence of the SR is one of the most 
intriguing problem in this area. Initially defined as a problem of "bona fide" stochastic 
resonance[4], it has created a vivid discussion in scientific literature[5, 4]. The main 
question in this discussion can be formulated as follows. It is a common knowledge, that 
the SR phenomenon manifests itself as a maximum of the signal-to-noise ratio (SNR) 
at the temperature, for which the Kramers rate equals the double driving frequency. 
However it is not possible to display a similar resonant characteristic as a function of 
frequency. 

On the other hand, experimental observation of the resonant escape from the 
potential well [6] has stimulated a number of theoretical papers [7, 8]. It has been 
shown that an underdamped particle driven by both noise and an AC force can perform 
a resonant escape from the potential well at a certain resonant frequency. However, the 
narrow resonant peak has been obtained only for underdamped or moderately damped 
systems, and only recently Pankratov and Salerno in [9] have shown that the mean 
transition time of an overdamped particle over a periodically oscillating potential well 
can exhibit a resonant dependence on the frequency of this oscillation in the case of 
strong (comparing to the depth of the well) driving forces. 

The problem we would like to present in this paper is whether an enhancement of 
the well-to- well transfer (jump) rate is possible in the overdamped bistable system for 
an arbitrarily strong (or weak) AC force. We believe that the answer to this question is 
also important because the unperturbed overdamped system does not have oscillatory 
dynamics, and, therefore does not possess oscillatory frequencies. We are also interested 
in finding the frequency dependence of the average jump rate for different values of noise 
and driving amplitude. 

2. The model 

Consider the overdamped particle in a double-well potential V(x) and under the action 
of a periodic external force and a white noise. The time evolution of the particle position 
x(t) is governed by the Langevin equation 

7 x = F(x, t) + £(t), F(x, t) = -V'{x) + E(t). (1) 

Here the dot stands for differentiation with respect to time, 7 is the dissipation 
(damping) parameter, and the double-well potential is chosen as in the 4 model: 
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V(x) = (x 2 — l) 2 /4. The relaxation time scale is given by t re i ax = 27 /V"(±l) = 7. 
The driving force E(t) = E(t + T) has been taken as 

E(t) = A tanh [// sm(ut)]/ tanh /x. (2) 

The parameter /x controls the shape of the function E(t) (see Fig. 1 for details). In the 
limit /i = we obtain simple sinusoidal drive, E(t) = Asin(ujt). In the opposite limit, 
/x — > 00, the driving force is of the stepwise form, e.g. E(t) = +A or —A during each of 
the halfperiods of the driving force. The importance of the shape of the driving force 
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Figure 1. Shape of the driving force E(t) for different values of the parameter /x: 
/i = fO (curve 1), yU = 2 (curve 2) and it = (curve 3). 

will be discussed later. The AC force is assumed to be weak, e.g. A <^ 1, so that the 
double-well system does not turn itself into a single-well system during the oscillations 
of E(t). 

The temperature effect is given by a zero-mean stochastic force (£(£)) = and 
auto-correlation function 

mm) = 2I>y5(t-lf), (3) 

with D being the strength of the noise (dimensionless temperature). Similarly to 
Eckman and Thomas [10] we are going to characterize the phenomenon by the average 
number of jumps (or the jump rate) from one well to the other one: 

n = lim — , (4) 

T^OO f 

where Nj is the number of jumps during the time interval r. 

In the absence of the periodic driving (A = 0) the average jump rate is defined by 
the well-known Kramers formula (see, for example, [11]). Note that for our choice of 
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the potential the height of the barrier is AV = V{x max ) — V(x min ) = 1/4, and, as a 
result, the Kramers rate equals 



We introduce the amplification factor 77 as a ratio of the average number of jumps in 
the presence of the driving fields to the jump rate when it is absent: i](uj) = n(uu)/rk- 
Since it is not possible to solve the corresponding Fokker-Planck equation analytically 
for all ranges of the driving frequency, we will tackle our problem numerically. 

3. Resonant nature of the jump rate 

We integrate the corresponding Langevin equation (1) with the Runge-Kutta method 
and compute numerically the number of particle jumps from one well to another one. We 
assume that the jumping event took place when the particle had reached the opposite 
well of the potential. The numerical experiment has been performed during the long 
times r ^> T = 2n/uj,l/r k , t re i ax = 7, so that the time averaging will eliminate any 
possible dependence on the initial phase. 

3.1. Main results of the numerical simulations 

Numerical simulations clearly show resonant behaviour of the amplification factor rj{uj). 
In Fig. 2 we show the dependence 77 on the driving frequency for the different values of 
A and D. Each of the four panels in the figure corresponds to the certain ratio A/D. 
The resonant nature of the phenomenon can be seen for all values of A/D, however it 
attains considerable values if A/D > 1. For example, in the panel (d) we present the 
case of A/D = 5 where the enhancement can reach the values of 10 -j- 15. Note that we 
did not continue down to the values 00 < 10~ 4 since it will require too long averaging 
times. However, even in the case (d) maximum of the dependence 77 = i](uj) can be seen. 

Particulars of the dynamics for different values of the driving frequency are shown 
in Fig. 3. The upper panel shows the time evolution of the particle in the case when 
no periodic driving force is present. Here jumping events are extremely rare and occur 
at the rate of order of r k ~ 8 x 1CT 7 [given by the Kramers formula (5)]. The middle 
panel shows the particle dynamics when the driving frequency equals u = 0.001, that 
corresponds to the frequencies close to the resonance. Gross enhancement of jumping 
events is clearly seen. The last panel deals with the case of frequencies much higher 
than the resonant frequency (cu = 0.5 in this case) where the intensity of jumping events 
is again low. 

3.2. Resonance in the jump rate as a switching between different time scales 

There are three important time scales in our system: the period of the driving force 
T = 1-r/lo, the mean residence time in the well, inverse to the Kramers rate, l/r k and 
the relaxation time of the system, t re i ax = 7 1/Vfc- It is not difficult to see that 




(5) 
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Figure 2. Dependence of amplification of the jump rate -q on the logarithm of the 
driving frequency for \i = 0, 7 = 1, and different ratios A/D. Panel (a): A/D = 1/2; 
D = 0.04(o),L> = 0.06(»). Panel (b): A/D = 1; D = 0.04(c), D = 0.06(»). Panel (c): 
A/D = 2; D = 0.03(»), ^ = °- 04 13 = °- 06 (+)• Pancl ( d ) = 5, £> = 0.02. 




Figure 3. Three examples of time evolution of the dynamic variable x(t) for 7=1, 
(j, = 0, D = 0.02 at A = 0, w = (upper panel); A = 0.1, w = 0.001 (middle panel) 
and w = 0.5 (lower panel). 



the dependencies 77 = 77(0;), shown in Fig. 2, consist of three distinct parts: (i) f](uj) 
increases in the adiabatic limit {u <^ r k ) from zero frequencies and continues to grow 
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when iq exceeds the value of the Kramers rate r k ; (ii) the plateau in the intermediate 
area r k < u < I/7, where the change of i] is slow and (iii) the high frequency limit, 
00 ^> I/7 where r] decreases with uj until rj reaches unity. We shall try to explain the 
observed resonant enhancement of the jump rate the switching of the above timescales. 

First, we try to explain why the jump rate grows in the adiabatic limit. Initially 
we consider the case of sinusoidal drive (// = 0). In the adiabatic limit n — > r k only for 
very small driving amplitudes (A — > 0). Momentary escape rates from the lower (r_) 
and upper wells (r + ) have been calculated in [12]: 



r±(<f>) ^ r fc £±(A,0)exp 



, Asin0^ 



3A 

Z±(A,<f>) = l^ — sin0, (6) 



and are, in principle, functions of the the phase = ujt. The momentary residence 
times t±(4>) = l/r±(0) obviously also depend on the phase. It is easy to see that for 
the same phase the increase of the average residence time in the lower well does 
not compensate the decrease on the average residence time in the upper well In 
other words, i_ + t + 7^ 2/r&, and this "asymmetry" in average residence times must 
depend on the driving amplitude A in the nonlinear fashion. Using Eq. (6), we show 
the following: 

2 _ 2r_(0)r + (0) _ r (1 - ggf^j exp(-M^) 

i / A sin A\\ jKAs\t\A\ • 1 / 4 ein \ " ^ ' 



t_(0) + t+{4>) r_(0) + r+(0) C0S h (^) - M|l^ sinh 

After some simple manipulations one can obtain that for any physically reasonable 
values A and D and for any except = 0, ±7T the expression multiplied by r k in 
Eq. (7) will be less then 1. Thus, since the inequality i-(0) + t+(4>) 2/ r k holds, 
the same should hold for the average values: t_ + t + 7^ 2/r^. The average number of 
jumps from well to well can be treated as the inverse mean residence time in the well, 
t — (t- + 1 + )/2 = l/n. Strictly speaking, to compute the mean residence times in the 
upper and lower wells, i±, one must perform proper averaging over the phase along 
the interval < < n. Notice that the greatest contribution to the asymmetry in 
the jump rates comes from the times when the driving force E(t) attains its extreme 
values, t = nT/4, n = 1,3,5,.... And, on contrary, at the zeroes of E(t) [t = nT/2, 
n = 0,1,2,...] the potential is symmetric and their contribution to the asymmetry of 
residence times is minimal. Thus, the following estimate of the average jump rate (or 
of the amplification factor 77) is possible: 



'1 - 



9A2 VW 3A2 



= 1>V> r , x t^t = 7 v — v / \ .(8) 



r k [t.(0)+t + {0)] r fc [Mf)+Mi)] cosh sinh (±) 

Thus we can conclude that r](uo — >• 0) < I. Also, it is straightforward to show that 
rj — > 1 if A — > and decreases with growth of A. This is in accordance with numerical 
results of Fig. 2, where in panels (a)-(c) it is clearly seen that rj < 1 when the driving 
frequency tends to zero. Also from the numerical results of Fig. 2 one notices that with 
the growth of the driving frequency the amplification ratio starts to increase. In the 
adiabatic limit (when T — > 00) the asymmetry of the residence times in the upper and 
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lower wells is maximal. As long as we increase the frequency (decrease T), less jumping 
events can take place during one driving period. As a result, more often the jumping 
events will take place around the zeroes of E(t), and, consequently, the contribution of 
the extreme values of E(t) into the asymmetry between residence times must decrease. 
This must cause the growth of the jump rate. However, the most remarkable fact is that 
the amplification ratio r\ does not stop growing when it reaches unity, but continues to 
increase. 

Another limit of our problem is the limit of high driving frequencies 00 ^> I/7. In 
this case, as opposing to the adiabatic lot of oscillations of the driving force can 

take place in between the two consecutive jumping events. Thus, the system does not 
"feel" the AC force on average, and, since (E(t)) t = 0, the influence of the AC drive 
diminishes. Finally, in the limit 00 — > 00 the jump rate must tend to the Kramers rate. 
This is in accordance with previous studies [12, 13] of the resonant escape from the 
potential well. 

Thus, we can conclude that there must exist an optimal driving regime for which 
the Brownian particle jumps between the wells with the maximal intensity, and this 
regime must be somewhere in between the adiabatic and high-frequency limits, e.g. in 
the interval r k < uo < I/7. We have estimated numerically the value of the resonant 
frequency as tu res ~ and this value is much larger then the Kramers rate. In 

Fig. 4 we investigate this dependence in more detail. In the panel (a) we observe that 
the maximum of the dependence 77 = rj(u) shifts to the right when 7 is decreased (while 
all other parameters are fixed). In the panel (b) we present the behaviour of the resonant 
frequency u res and we see that it does not show any clear dependence on the amplitude. 
Although we do not have a rigorous theory to support such a formula, qualitatively it 
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Figure 4. Panel (a). Dependence of the scaled amplification factor on the driving 
frequency at A — 0.06, D = 0.03 and n = for different values of damping: 7 = 1 (•), 
7 = 0.5(o), scaled by 90%, and 7 = 0.1 (squares), scaled by 80%. 
Panel (b). The resonant frequency a function of (r k hfl 2 for A = 0.02 (o), 

A = 0.03 (+), A = 0.04 (*), A = 0.06 (•), A = 0.08 (o) and A = 0.1 (squares). The 
dashed line is the bisectrix of the straight angle and is shown as a guide to an eye. 
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is easy to understand that this behaviour of the resonant frequency is a result of two 
competing effects: increasing of the temperature tends to activate the jumping process 
while damping works against it. 

From this estimate of u res we see that the resonant frequency is linked with the 
value of the Kramers rate. Doubling the temperature may increase u res by an order of 
magnitude. Moreover, the steep growth of the curve w (cu) depends solely on the Kramers 
rate, as it is very clearly demonstrated in Fig. 2c for three different values of D and in 
Fig. 4a for three different values of 7. The position of the inflection point in the curve 
rj(uj), which approximately separates the areas of the steep growth and the plateau, 
moves to the right when the temperature D is increased. This inflection point is always 
greater then but is approximately the same order of magnitude as rj,. On the other 
hand, decrease of the dependence i](uj) is controlled only by the damping parameter. It is 
easy to see that the inflection point which separates the plateau from the high-frequency 
area does not depend on the temperature. However, it depends strongly on damping 
7 (see Fig. 4) and matches approximately the relaxation frequency l/t re i ax = 1/7 (see 
Figs. 2). This can be explained in the following way: the most probable dynamics is 
that the particle jumps from the upper well to the lower well. If the time during which 
the particle moves from one well to another one is of the order of T/2, the particle will 
appear again in the upper well and will perform a jump with a higher probability again. 
When u increases, such synchronization with the relaxation frequency will be lost and 
if to I/7 the particle will have no time to relax in between the oscillations of the 
potential. Thus, the jumping events will be more rare, the jump rate will decrease. As 
we enter the high-frequency area the system will react on the external periodic field very 
weakly. We can conclude that in order to get the maximal jump rate, the frequency of 
the drive must be synchronized with both the Kramers rate r k and the characteristic 
relaxation frequency I/7. 

3. 3. Resonant jump rate as a nonlinear response 

An important question is how the observed phenomenon depends on the amplitude of 
the driving force, A. In this subsection we show that this dependence has the nature of 
nonlinear response to the driving force. In principle, there are many ways to demonstrate 
this. As one example, in Fig. 5 we have plotted the dependence of the amplification 
factor at its maximal value as a function of the ratio A/D with D being fixed for each 
of two curves, shown in the figure. Both curves have been fitted with the parabolic 
functions. For the case D = 0.04 we have obtained 1 + 0.245(A/ D) 2 and for D = 0.06 
we have obtained l + 0.21(A/D) 2 . Another example which supports the above statement 
about the nonlinear response is the behaviour of the lower boundary of the amplification 
ratio rj in the adiabatic limit [see Eq. (8)]. The expansion of this formula in the powers 
of A/D yields 1 - (17 / 16) (A /Df + O ((A/D) 4 ). These results unambiguously show that 
resonant enhancement of the jump rate can not be described within the framework of 
the linear response. 
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Figure 5. Dependence of the amplification ratio at the resonant value of the driving 
frequency on the ratio A/D for 7 = 1 and fi = 0: D = 0.04 (*) and D = 0.06 (o). 
Solid lines are the best fits by the parabolic function (see text for details). 

3.4- Dependence on the shape of the driving field E(t) 

Finally we discuss the generality of the obtained results with respect to the shape of the 
driving force E(t) since the distribution of the residence times t± on the initial phase has 
been omitted when the adiabatic limit has been considered. To check this, numerical 
simulations have been performed for the two cases of the driving field E(t). In the 
Fig. 6 we show the dependence of the amplification ratio when the drive is sinusoidal 
(/i = 0) and when it is almost stepwise (/z = 10). In the latter case one can assume the 
distribution of the residence times along the halfperiod of the oscillation of the driving 
force to be independent on the phase 0. No qualitative differences can be observed for 
these cases, only the quantitative ones. In the adiabatic limit, the amplification factor 
is less then one and increases with the growth of the driving frequency regardless of the 
shape of the driving field E(t). 

Another interesting observation is that the maximal amplification of the jump rate 
depends on the shape of the driving field. We see that for the stepwise shape (large 
fi) we obtain much better amplification. If E(t) is of the stepwise shape, the time 
intervals when the potential wells are almost not desymmetrized, are very short. In 
other words, almost all the time one well is in the uppermost position while another 
one is in the lowest. The time intervals when the potential is almost symmetric do not 
contribute anything into the observed resonance, since in that case the jump processes 
occur with the rate close to r^. Thus minimization of these time intervals effectively 
must enhance the jumps and selecting appropriately the shape of E(t) one might increase 
the amplification significantly. However, this topic is beyond the scope of our paper. 

4. Conclusions 

To summarize our results, we have found that the jump rate of an overdamped Brownian 
particle in the double-well potential, driven by a periodic force with zero mean depends 



Resonant enhancement of the jump rate in a double-well potential 



10 




_4 _3 -2-1 1 

log CO 



Figure 6. Dependence of the amplification of the jump rate rj on the logarithm of the 
driving frequency for 7 = 1, D = 0.03 and A — 0.09 for different shapes of the driving 
field E(t): fi = (o) and ^ = 10 (+). 

resonantly on the frequency of the drive. This phenomenon exists for any ratio of the 
driving amplitude A to the noise strength, D, however is very well pronounced if A > D. 
We also have observed that this phenomenon occurs for different shapes of the driving 
force E(t). For the stepwise function E(t) the jump rate can be twice larger then in 
the case of the sinusoidal drive with the same amplitude. An interesting result has been 
observed in Ref. [14], where diffusion of an overdamped particle in the spatially periodic 
potential has been studied. The particle is also under the influence of the time-periodic 
force. Among other results, authors have noticed "acceleration" of diffusion for certain 
values of the driving period. However, while it was only in the limit A/D ^> 1 they 
observe the effect, in our case it happens for any ratio of A/D. 

We also would like to mention some analogy of the observed phenomenon with the 
resonant activation over a fluctuating barrier, which can be driven either randomly as in 
Ref. [15] or deterministically via periodic force (see Ref. [9]). These papers have shown 
that the mean first-passage time (MFPT) as a function of the intensity (rate, frequency) 
of the barrier fluctuations displays a local minimum. In our problem the driving force 
E{t) effectively induces the oscillating change of the depth of the wells and the average 
number of jumps n is nothing but an inverse MFPT. 

Finally we would like to stress the importance of the overdamped case since the 
dynamics of many biological objects is modeled by overdamped systems, so the above 
problem could have a wide range of applications. 
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